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ABSTRACT 

We explore the near-infrared properties of galaxies within 27 galaxy clusters using data from the 
Two Micron All Sky Survey (2MASS). For a subsample of 13 clusters with available X-ray imaging 
data, we examine both the properties of the galaxies and the intracluster medium. We show that 
the i-T-band luminosity is correlated with cluster mass, providing a binding mass estimate accurate to 
45%. The mass-to-light ratio in our ensemble increases by a factor of ~ 2 over the cluster mass range 
(1O 14 M0 — 10 15 Mq). We examine the total baryon fraction, showing that it is an increasing function of 
cluster mass. Using the mass-to-light ratio of massive clusters, we find that SIm = 0.19 ±0.03; using the 
total baryon fraction we find that — 0.28 ± 0.03, in good agreement with recent cosmic microwave 
background (CMB) anisotropy constraints. Differences between these two estimates suggest that the 
-ftf-band mass-to-light ratio in massive clusters may be lower than that in the universe by as much as 
-30%. 

We examine the stellar mass fraction, the intracluster medium (ICM) mass to stellar mass ratio and 
the cluster iron mass fraction. The stellar mass fraction decreases by a factor of 1.8 from low to high 
mass clusters, and the ICM to stellar mass ratio increases from 5.9 to 10.4 over the same mass range. 
Together, these measurements suggest a decrease of star formation efficiency with increasing cluster 
mass and provide constraints on models of the thermodynamic history of the intracluster medium. The 
cluster iron mass to total mass ratio is constant and high, suggesting that some efficient, and uniform 
enrichment process may have taken place before the bulk of stars in cluster galaxies formed. 

Subject headings: cosmology: observation - galaxies: clusters, luminosity function - infrared: galaxies 



1. INTRODUCTION 

Galaxy clusters are central to many cosmological stud- 
ies. The local cluster abundance is a direct measure of 
the present-day matter power spectrum amplitude and the 
dark matter density (White et al. 1993a; Viana & Liddle 
1999; Reiprich & Bohringer 2002). The cluster redshift dis- 
tribution constrains the growth of density perturbations, 
and this effect can be used to constrain the amount and 
nature of the dark matter and dark energy (Wang & Stein- 
hardt 1998; Haiman et al. 2001; Holder et al. 2001; Levine 
et al. 2002; Majumdar & Mohr 2002). The baryon fraction 
of the most massive local galaxy clusters has been used 
to constrain the matter density parameter (White et al. 
1993b; David et al. 1995; White & Fabian 1995; Mohr et al. 
1999; Grego et al. 2001), and the spatial distribution of 
galaxy clusters has been used to measure the power spec- 
trum of density perturbations (Bahcall & Soneira 1983; 
Miller & Batuski 2001; Schuecker et al. 2002). 

The galaxy cluster halo mass is crucial to the interpreta- 
tion of each of the cosmological measurements mentioned 
above. The cluster halo mass, which is dominated by dark 
matter, is not a direct observable; a good mass indicator is 
required to bridge the gap between the observational and 
theoretical realms. Detailed mass estimators relying on 
hydrostatic or virial equilibrium or on the weak lensing dis- 
tortions of background galaxies have been developed, but 



generally these estimates require extensive X-ray and op- 
tical imaging and/or spectroscopy. The large cluster sur- 
veys that are being planned and carried out to measure the 
cluster redshift and spatial distribution will produce large 
samples dominated by systems that are near the detection 
threshold. Deep follow up of each of these systems with 
an X-ray telescope or multiobject spectrograph is implau- 
sible, and so the analyses of these future surveys will rely 
on the existence of simpler and perhaps less accurate mass 
estimators. Specifically, these surveys will employ scaling 
relations between simple cluster observables like the X-ray 
luminosity or the Sunyaev-Zel'dovich effect (SZE; Sunyaev 
& Zel'dovich 1970, 1972) luminosity and cluster halo mass. 

Observations over the past several years have demon- 
strated that clusters exhibit significant regularity (Mohr 
& Evrard 1997; Arnaud & Evrard 1999; Mohr et al. 
1999), and we now have several locally calibrated mass- 
observable relations. Among these are the mass-X-ray 
temperature relation (M — Tx, e.g. Finoguenov et al. 
2001b), the mass-luminosity relation (M — Lx, e.g. 
Reiprich & Bohringer 2002) and the mass-velocity disper- 
sion relation (M — ay, e.g. Girardi et al. 2000). Here we 
use the Two Micron All Sky Survey (2MASS) to explore 
whether the galaxy light itself is a good indicator of the 
cluster mass by examining the cluster halo mass-if-band 
light relation. Such a relation would be extremely valuable 
for interpreting any cluster survey with ncar-IR (hereafter 
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NIR) coverage. 

In addition to the cosmological studies, there is a range 
of unresolved questions surrounding the star formation his- 
tory of the universe, and the thermodynamic history of 
the intergalactic and intracluster medium (hereafter ICM). 
The galaxy NIR light in combination with X-ray imaging 
and spectroscopy are well suited to address these issues, 
because the X-ray observations provide measures of the 
halo and ICM masses and the NIR light appears to be 
a good tracer of the total stellar mass. The effects of 
dust in NIR wavebands are rather small compared to op- 
tical and ultraviolet (UV) wavebands. The NIR light is 
less sensitive to the recent star formation history in galax- 
ies. Studies of a large sample of spiral galaxies (Gavazzi 
et al. 1996) indicate that the iJ-band luminosity is pro- 
portional to the galaxy dynamical mass, implying that H- 
band light can serve as a better mass estimator than U, 
B or V^-band light. Modeling suggests that the isT-band 
mass-to-light ratio can still vary by as much as a factor 
of two over a range of galaxy Hubble type, color, and star 
formation histories (Madau et al. 1998; Pahre et al. 1998), 
so care must be taken to account for these variations wher- 
ever possible. 

Massive galaxy clusters largely exhibit self-similar scal- 
ing relations, but that behavior does not hold for less mas- 
sive clusters and groups (e.g. Lloyd-Davies et al. 2000; Voit 
et al. 2002). Non-gravitational processes, such as feedback 
due to star formation activity in galaxies or AGN heat- 
ing may significantly alter the structure of the lower mass 
systems, but important questions remain about which pro- 
cesses are important. By studying the state of the baryons 
in present epoch clusters, one can in principle unravel de- 
tails of their formation and evolution. For example, the 
distribution of the ICM gas may give clues to the impor- 
tance of non-gravitational heating (Cavaliere et al. 1997; 
Mohr et al. 1999; Bialek et al. 2001); the total stellar con- 
tent within a cluster measures the cooled baryon fraction, 
which in turn describes the efficiency of galaxy formation 
in the clusters (Bryan 2000); the metallicity of the ICM gas 
also provides stringent constraints on the star formation 
history of cluster galaxies, as well as on the mechanisms 
that transport the metals from the galaxies into the ICM 
(Arnaud et al. 1992; Renzini 1997). Among the 27 clusters 
in our sample, we have measured the ICM mass for a sub- 
sample of 13 systems. This, together with the total stellar 
mass inferred from the NIR luminosity, gives us an oppor- 
tunity to address different preheating, star formation and 
enrichment scenarios. 

This paper is structured as follows. In §2 we describe 
how we infer the cluster virial mass from X-ray tempera- 
ture, how the total NIR galaxy luminosity is obtained from 
the 2MASS data, and how the ICM mass is measured us- 
ing X-ray data. We present the NIR mass-luminosity re- 
lation, as well as the mass-to-light ratio for our sample of 
clusters in §3. In the rest of §3 we examine constraints on 
the density parameter using the cluster mass-to-light 
ratio (§3.1) and the total baryon fraction (§3.2). We com- 
pare the results and discuss the cosmological implications 
in §3.3. In §4 we examine the cluster baryon reservoirs 
in more detail. The discussion includes the star forma- 
tion efficiency on cluster mass scales (§4.1), the relative 
distribution of baryons in the ICM and in galaxies (§4.2), 



and the ICM iron abundance and clues on the enrichment 
(§4.3). Finally, we present our conclusions, discuss the 
robustness of our analysis, the possible systematics, and 
prospects for further investigation in §5 & §6. 

Throughout the paper we assume the density parame- 
ters for the matter and the cosmological constant to be 
VLm = 0.3, Oa = 0.7, respectively, and the Hubble param- 
eter to be H = 70h 7a km s _1 Mpc -1 . 

2. JOINT NIR AND X-RAY ANALYSIS 

We use the observed X-ray mass-temperature relation 
together with published X-ray emission weighted mean 
temperatures, 2MASS second incremental release NIR 
data, and X-ray imaging data to study trends in the NIR 
and X-ray properties of galaxy clusters. Our cluster sam- 
ple is based on existing cluster samples (David et al. 1993; 
Mohr et al. 1999, hereafter MME; Reiprich & Bohringer 
2002) and is limited by the incomplete sky coverage in the 
2MASS second release (~ 47% of the sky). We return to 
the issue of sample selection and possible effects on our re- 
sults in §2.2.2 & 5. The relevant cluster parameters (X-ray 
emission center, redshift z, X-ray emission weighted tem- 
perature Tx) are gathered from the above references and 
others (Edge & Stewart 1991; White 2000). We only con- 
sider the clusters that have reliable Tx measurements and 
are reasonably beyond the galactic plane (|6| > 10°). The 
2MASS second release contains 27 clusters for which lumi- 
nosity functions can be estimated in the method described 
below. Among these 27 clusters, the ICM masses Micm 
of 13 were measured by MME, based on an analysis of 
ROSAT PSPC observations. We shall refer to these as 
the MME subsamplc. A list of our cluster sample ap- 
pears in Table 1. In short, our cluster sample is at low 
redshift (0.016 < z < 0.09 with mean redshift - 0.05), 
and spans a range of X-ray emission weighted temperature 
(2.1 < Tx < 9.1) that corresponds to about an order of 
magnitude in cluster binding mass (0.8-9xl0 14 /if 1 M Q ). 

2.1. X-ray Estimate of Cluster Mass 

We use an observed X-ray M — Tx relation to calcu- 
late the binding mass for each cluster. Finoguenov et al. 
(2001b) provide several M — Tx relations that arise from 
different subsets of a cluster ensemble. We use the relation 
obtained by fitting to the clusters with ASCA temperature 
profiles that are more massive than 3.57 x 10 13 h^Q M Q ; all 
our clusters lie in this mass range. To be definite, we use 

1 co + O-06 

M / TV \ ° -0.07 

where M 50 o is the mass enclosed by r 50 o , within which the 
mean overdensity is 500 times of the critical density of the 
universe p c . Because our clusters are all low redshift, we 
neglect the effects of cosmic density evolution. 

Within our fiducial cosmological model (Qm = 0.3 and 
Oa = 0.7) we calculate the cluster radius r 50 o and its an- 
gular extent #500 = r^oo/dA, where cIa is the angular diam- 
eter distance. We analyze the galaxy and ICM properties 
within this region of the cluster. At the low redshift range 
spanned by our cluster sample the sensitivity to cosmology 
in the angular diameter distance is unimportant. 

2.2. NIR Luminosity from 2MASS Data 
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Knowing the cluster angular extent, we search the 
2MASS extended source catalog to identify cluster mem- 
ber galaxies. Among the three bands (J, H and K 8 ) 
2MASS provides, K s has been studied most extensively. 
For this reason we shall only consider the AT s -band (here- 
after denoted as AT-band for simplicity) in this paper. We 
use the "default" K s magnitudes. The official 2MASS sec- 
ond release completeness limit for the default AT s -band 
is Ku m = 13.5 5 . However, after visual inspection of 
the magnitude distribution of the galaxies, we choose 
Kum = 13.3 as our completeness limit. 



II II II 


iii 

I 

T } 


III 

{ t ' 


iii 

< 


1 

i 

! 


i 




i 

iii 


i 


5 


10 


15 


20 


25 



Cluster § 

Fig. 1. — A comparison between the background galaxy surface 
density as obtained using the statistical method (H a tat, solid line), 
and the "annulus" method (Ej, 9n , solid points, see text) 

2.2.1. Background Correction 

To correct both the number counts and fluxes from back- 
ground (and foreground; hereafter we refer to all non- 
cluster galaxies as background galaxies) galaxies in esti- 
mating the cluster NIR luminosity, we use the observed K- 
band galaxy brightness distribution d Y> sta t/dK (Kochanek 
et al. 2001). Integrating this quantity with respect to mag- 
nitude to our adopted 2MASS completeness limit gives the 
statistical galaxy surface density £ s tat- This, multiplied by 
the cluster angular area 7T0| OO , gives the expected number 
of background galaxies. The expected flux from the back- 
ground galaxies is F s t a t — ^foo 

(dZ stat /dK)FdK, 

where F is the galaxy flux. 

As a check of this statistical method, we also perform 
an "annulus" background correction; that is, for an as- 
sumed galaxy distribution profile with unknown normal- 
ization and a unknown constant background surface den- 
sity, for each cluster we measure the number of galaxies 
within 05oo (N cen ) and within an annulus outside the clus- 
ter virial region (N ann ), whose inner and outer radii are 
1.5 0200 and 3 9 2 oo, respectively, where 8 2 oo = ^oo/cU ■ We 
convert between r2oo and r5oo using the "universal" dark 
matter distribution profile found in N-body simulations 
(Navarro et al. 1997, NFW); we use a fixed "concentra- 
tion" parameter c = 5, and in that case r2oo = 1-51 rsoo- 

5 http:/ /www. ipac.caltech.edu/2mass/releases/second/doc/explsup. 



Background and cluster galaxies lie in both the central 
and annular regions. Using the observables N cen and N ann 
and the NFW model for the distribution of galaxies, we 
estimate the normalizations of the cluster galaxy distribu- 
tion and the uniform background galaxy surface density 
Sfc grl . The number of cluster galaxies within $500 is then 
A^cen ^bgn ^ ^loo- For some clusters the annulus is only 
partially surveyed in the 2MASS second release, and we 
correct for this partial coverage. In some cases there are 
obvious galaxy clusterings (from visual inspection) in the 
annulus region; we avoid this possible contamination of 
the background annulus by other clusters or groups by ex- 
cluding the area occupied by the clustering in the annulus. 

Fig 1 contains a comparison of ~S s tat determined from 
the statistical approach and £{, gn as determined from the 
"annulus" method. There is general agreements between 
£ s tat and Sf, 9re , indicating that local and statistical back- 
ground corrections lead to the same results. For the analy- 
ses that follow, we adopt the statistical background surface 
density in making background corrections to the cluster 
galaxies. Note that, for 3 clusters, determined from 
the "annulus" method is negative, which may either be due 
to the correction to the partially surveyed annulus region, 
or be due to a significant mismatch between the actual 
galaxy distribution in these clusters and the NFW model 
we are using. We have checked all the scaling relations pre- 
sented in the following analysis using the local background 
estimate, and the changes to slope and normalization all 
lie within the la uncertainties. 
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Fig. 2. — Measured M* and </>* for our cluster sample (see Table 
1). Solid circles are clusters with kT > 4keV. Larger symbols de- 
note clusters with z < 0.05, while smaller symbols refer to z > 0.05 
clusters. 



2.2.2. Luminosity Function 

We transform the apparent magnitudes into solar lu- 
minosities using Kq = 3.39 (Kochanek et al. 2001), the 
cluster distance and Galactic extinction: 

M K = K - 25 - 5 log 10 (dL/lMpc) - A K , (2) 
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Table 1 

Basic Descriptions of the Clusters 



(1) 


(2) (3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 








M, - 5 log h 70 


-^500 


T500 


/star 


flCM 


Name 


z (koV) 


W S ai (/4 MpC- 3 ) 
r ^1.81 


(mag) 


(/^ 2 IO^Lq) 


(h 70 T Q ) 

^r!4 


(io- 2 ) 

-, ,tU.4H 


(IO" 2 ) 



A119 

A262 

A496 

A644 

A754 

A1367 

A1651 

A2255 

A2319 

A3266 

A3558 

A4038 

A133 

A548e 

A1185 

A1650 

A1767 

A2107 

A2147 

A2151 

A2440 

A2589 

A2593 

A2626 

A2634 

A3389 



0.0440 
0.0161 
0.0328 
0.0704 
0.0528 
0.0216 
0.0850 
0.0800 
0.0564 
0.0594 
0.0480 
0.0283 
0.0569 
0.0410 
0.0304 
0.0845 
0.0701 
0.0421 
0.0351 
0.0370 
0.0904 
0.0416 
0.0433 
0.0573 
0.0312 
0.0265 



5.8" 



1.36 

°- yl 0.04 
6.59«; 

8 5 03 

6 3°' :3 ^ 
6.87"| 

9 12 

6 2 oS 

iSsg 3 

M 6 
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3.8, 
3.1! 
3.9' 

4.31°;* 

2.4f 
3.8.' 
3-7,,. 
3 l 1 5 
2^ :g 

o 7 0.18 
°- '0.18 
2 1 0S> 

z - 1 n.6 



o.n'o 
0.06 
0.16 
0.14 



39 
41 
40 
17 
44 
59 

8 
16 
61 
30 
61 
32 
13 
31 
27 

6 

7 
22 
42 
23 

7 
19 
23 
14 
27 
27 



3.80| 
7.312 
7.41 
5.57? 
4.04* 
7-25 2 24 
17.55*$ 
3.36, §fi 
5. 95j 
4.92* 
7-59 
6.50^-; 

7.71? 
4.01? 
4.66 2 

3.61, 
5.69f 
5.32} 
3.11} 
6.74| 
7-17f. . 

7 58 1156 

' - 00 5.46 
5 IS 1 - 83 

°i-§a 



.is; 

8.275 



-25.24 ± 0.20 
-24.69 ± 0.27 
-24.59 ± 0.20 
-24.68 ± 0.19 
-24.86 ± 0.16 
-24.77 ± 0.20 
-24.06 ± 0.17 
-25.30 ± 0.21 
-24.83 ± 0.12 
-25.03 ± 0.18 
-25.09 ± 0.14 
-24.68 ± 0.24 
-23.95 ± 0.22 
-25.15 ± 0.22 
-24.81 ± 0.37 
-24.55 ± 0.26 
-24.84 ± 0.39 
-25.04 ± 0.28 
-24.62 ± 0.19 
-25.87 ± 0.29 
-25.79 ± 0.34 
-24.29 ± 0.29 
-24.91 ± 0.28 
-24.95 ± 0.34 
-24.53 ± 0.26 
-25.05 ± 0.31 



8.32 ± 0.05 
2.51 ± 0.04 

4.96 ± 0.03 
8.45 ± 0.03 
10.99 ± 0.03 
4.42 ± 0.03 
13.89 ± 0.01 
9.81 ± 0.06 
17.30 ± 0.02 
9.94 ± 0.04 
13.91 ± 0.03 
3.23 ± 0.04 

4.97 ± 0.02 
5.72 ± 0.05 
3.25 ± 0.06 
6.74 ± 0.04 
3.96 ± 0.07 
4.53 ± 0.06 
5.19 ± 0.03 
5.18 ± 0.14 
6.15 ± 0.15 
3.39 ± 0.03 
4.56 ± 0.05 
4.70 ± 0.07 

3.33 ± 0.03 
3.29 ± 0.07 
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40 

44 

59 

68 



33 

54 
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45 
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67" 
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35^ 
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1.58" 

1.76° 

1.71 

1.33° 

1.17 

1.79° 

2.33°, 

1.45 

1.65° 

1.71 

2.72° 

1.53° 

1.79* 

2. 78°, 

1.12 

1.03° 

1.27* 

1.35° 

1.27° 

3.65} 

2.14° 

1.27° 

2.22} 

2.52 3 

1.24° 

2.82? 



12.03 ± 0.75 
10.17 ± 0.39 
13.47 ± 0.32 

13.04 ± 0.47 
13.90 ± 1.50 
11.85 ± 0.45 
13.69 ± 0.70 
13.14 ± 2.60 
16.69 ± 0.52 
16.66 ± 1.31 
15.83 ± 0.38 
11.64 ± 0.73 



Note. — Columes: (1) Name; (2) Rcdshift; (3) Emission— weighted mean temperature; (4) Estimated number of member galaxies; (5) Char- 
acteristic number density; (6) Characteristic magnitude; (7) Total luminosity within rsoo; (8) X-band mass-to— light ratio; (9) Stellar mass 
fraction; (10) ICM mass fraction. Uncertainties in all columns quoted at 1<t level. M* & </>* calculated assuming a = —1.1. All calculations 
done with h^o = 1. 



where <1l is the luminosity distance, and Ak accounts for 
the Galactic extinction. We ignore k-corrections, because 
our clusters are all at very low redshifts. The Galactic ex- 
tinction, which has a small effect on the overall luminosity, 
is assumed to be the same for the entire cluster; we obtain 
the extinction values from the NASA/IPAC Extragalactic 
Database . 

Using the NFW model with concentration c = 5, we 
deproject to convert the projected number and luminos- 
ity into an estimate of the number N a bs and luminosity 
L b s of galaxies within t^qq. The deprojection essentially 
reduces the background corrected luminosity and number 
by 20 — 30%. These quantities divided by the cluster vol- 
ume Vcis = (47r/3)r5 00 are the cluster mean number den- 
sity n b s and luminosity density j oos , which are related to 
the cluster luminosity function <j){L): 

/•oo />oo 

n obs = / 4>(L)dL and j obs = / Lcj)(L)dL, (3) 

JL lim JL Hm 

where Lu m is the luminosity corresponding to Ku m for 
each particular cluster. We use the Schcchter (1976) form 
for the luminosity function: 

(j) (L)dL = ^^y e- L / L ^d^y (4) 

where and <p* are characteristic luminosity and num- 
ber densities, respectively, and a is the faint-end power 
law index. 

In practice, the 2MASS data do not go deep enough to 
provide useful constraints on the faint end slope a for all 
of our clusters. Thus, we use external constraints to set 



a. The NIR cluster Schechter function faint end slope is 
generally found to be —1.3 < a < —0.8 (e.g. Mobasher & 
Trentham 1998; Andreon & Pello 2000; Tustin et al. 2001; 
Balogh et al. 2001a; De Propris et al. 2002). In the follow- 
ing we adopt a = —1.1. (We discuss the effect of different 
a on our results in §5.) Given n i, s , jobs, ct and Lu m , we 
solve for 0» and L* using Eqn 3. 

Fig 2 contains the derived 0» and M* (the absolute 
magnitude corresponding to L»). The distribution of M* 
(M» - 5\ogh 7Q = -24.88 ± 0.02) in Fig 2 and Table 1 
shows no clear correlation with cluster mass or redshifts, 
but many nearby clusters exhibit larger 0*. Our cluster 
sample is based on existing X-ray cluster catalogs, and 
so there is a tendency for the low temperature systems to 
lie at low redshift. For example, the mean redshift of the 
lower temperature half of the sample (i.e. kTx < 4 keV) 
is 0.040 ±0.001 as compared to 0.060 ±0.001 for the other 
half. This could potentially pose a problem, because any 
rcdshift related systematic could masquerade as a mass 
related trend in our cluster population. We probe for sys- 
tcmatics by examining the luminosity function parameters 
within the hot and cold cluster subsamples. Dividing each 
subsample at redshift 0.05, we find no statistically signifi- 
cant differences with redshift in M» or 0* for either of the 
subsamples. However, within the same redshift range, 0* 
is significantly higher for the low temperature subsample, 
consistent with the shallow light-mass relation presented 
above. 

In obtaining the Schechter parameters by using the ob- 
served light and number of galaxies, we do not consider the 
brightest cluster galaxy (BCG), which typically contains a 
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significant fraction of the cluster light and would bias the 
parameters of the cluster luminosity function. We account 
for the BCG light separately, modeling cluster luminosity 
functions with Schechter functions plus the light of the 
BCG. There are two clusters whose Schechter parameters 
appear to be different from that of other clusters: A133 & 
A1651. Their are relatively dim, and their about 
twice as large as others. This may be because of their 
relatively small number of member galaxies, and their ex- 
traordinarily bright central dominant galaxies: removal of 
the luminosity contribution from the BCGs dramatically 
reduces and boosts These two clusters underscore 
the fact that our integral approach to solving for and 
4>» leads to correlated errors in L* and that maintain 
an accurate estimate of the total cluster light even when 
estimates of and individually can be quite inaccu- 
rate. 



2.2.3. Estimating Total Light and Stellar Mass 

With the Schechter parameters and £*, we estimate 
the total cluster luminosity by integrating over the lumi- 
nosity function: L to t — Vcis L* 0* T(a + 2) + Lb cg . To 
reduce the dependence of our light estimate on the faint 
end slope parameter a we truncate the luminosity function 
at the absolute magnitude Mi ow — —20: 



j'tot - Lbcg/VcU + 



L<j>(L) dL. 



(5) 



For our sample and choice of a, j' tot differs from j to t only 
at 2 — 3% level. In the analyses that follow we make no 
distinction between these two, and use the "truncated" 
luminosity as our luminosity estimate throughout. 

We also apply a correction for the fact that the the de- 
fault 2MASS magnitudes underestimate the total light of 
the galaxies. By comparing the 2MASS default (isophotal) 
magnitudes of their galaxy sample with the total magni- 
tudes of deep photometry (Pahre 1999), Kochanek et al. 
(2001) found a systematic difference of 20%. We multiply 
Ltot by the same factor of 1.2. 

To obtain the total stellar mass, we multiply the total 
luminosity L to t by the average stellar mass-to-light ratio 
T s t ar (Tx) for each cluster. We estimate the mean stel- 
lar mass-to-light ratio from the observations of ellipti- 
cal and spiral galaxies (Gerhard et al. 2001; Bell & de 
Jong 2001), and take into account the varying spiral frac- 
tion as a function of cluster temperature (Bahcall 1977b; 
Dressier 1980, notice that the spiral fractions were deter- 
mined using £?-band light rather than A-band; we ne- 
glect possible wavelength dependence of the spiral frac- 
tion here). The details of our approach are presented in 
the Appendix. We show in Fig 3 the resulting mean stellar 
mass-to-light ratio T s t a r, and compare it with the result 
from the 2dF galaxy redshift survey (Cole et al. 2001), 
which was obtained by using the observed galaxy colors 
and redshifts to constrain their evolution, based on the re- 
cent version of a population synthesis code (Bruzual A. & 
Chariot 1993). This value 0.73T Q (where T Q is the solar 
A-band mass- to- light ratio) is consistent with our lowest 
temperature clusters (highest spiral fraction), but over the 
full range of cluster masses examined in this sample, the 
stellar mass-to-light ratio is expected to vary by ~10%. 
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Fig. 3. — The A"-band mean stellar mass— to-light ratio T s tar(Tx) 
based on the observed mass— to-light ratio in elliptical and spiral 
galaxies and spiral fractions in clusters. See Appendix for more de- 
tails. The dashed-line is the mean value from the 2dF survey (Cole 
et al. 2001) *£ s tar,2dF = 0.73 Tq, assuming a Kennicutt (1983) in- 
tial mass function. 



2.3. Estimating the ICM Mass from X-ray Emission 

We calculate the ICM mass for the 13 clusters in the 
MME subsample using fits to the ROSAT Position Sensi- 
tive Proportional Counter (PSPC) X-ray surface bright- 
ness fits from the literature (Mohr et al. 1999). The 
0.5:2.0 keV X-ray emissivity is relatively insensitive to 
temperatures (varying by 10% for a factor of two in- 
crease in temperature for clusters with temperatures above 
1.5 keV), so one needs not know the cluster temperature 
profile to determine the ICM density profile from an X-ray 
image. We have adjusted the density profiles to reflect im- 
proved temperature measurements in two of the clusters. 
Analysis of mock PSPC images of simulated galaxy clus- 
ters indicates that cluster ICM masses can be estimated 
with an accuracy of ~10% (Mohr et al. 1999) and that 
clumping in the ICM enhances the X-ray emissivity lead- 
ing to an estimated ~ 10% bias in ICM mass estimates 
(Mohr et al. 1999; Mathiesen et al. 1999). We correct for 
this bias by reducing ICM mass estimates for all 13 clus- 
ters by 10%. However, this reduction is counterbalanced 
when one considers the ~ 10% diminution of the cluster 
baryon fraction at rsoo compared to the global baryon frac- 
tion fifj/rjjvf as seen in numerical studies, where £lb is the 
total baryon density in the universe. We recalculate the 
ICM masses for this analysis using the M500 — Tx rela- 
tion adopted for the rest of the analysis (Finoguenov et al. 
2001b, see Eqn 1). 

3. NIR LUMINOSITY AND COSMOLOGY 

One critical issue for ongoing NIR cluster surveys is de- 
termining the degree to which the NIR light from the stel- 
lar populations within cluster galaxies traces the galaxy 
cluster binding mass. Cluster binding mass or virial mass 
is central to using the abundance of galaxy clusters to con- 
strain cosmology (although deep surveys contain enough 
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information to solve for cluster binding mass and cosmol- 
ogy simultaneously; Majumdar & Mohr 2002). If NIR 
galaxy light is a good indicator of cluster binding mass, 
then large solid angle surveys like 2MASS could be used 
to (1) define very large samples of clusters in the nearby 
universe (see Kochanek et al. 2002), and (2) provide direct 
constraints on a combination of the local power spectrum 
of density fluctuations and the mean matter density. Due 
to the limited size of our sample, we simply focus on the 
relation between the virial mass (inferred from X-ray ob- 
servations) and stellar luminosity (obtained from 2MASS 
photometry) of nearby clusters, discussing the implications 
for future NIR cluster surveys. In addition, we use the NIR 
light and measurements of cluster ICM masses to address 
the baryon fraction, the cluster mass-to-light ratio, and 
the cosmological matter density parameter. 

kT x (keV) 
2 4 6 8 10 
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Fig. 4. — The cluster K-band luminosity-mass correlation. The 
uncertainties in L500 are rather small, typically less than 10% (see 
Table 1). The solid line shows the best-fit relation of slope of 0.69 
(Eqn 6), while the dashed line shows a slope of 1.0 behavior. 

Figure 4 contains a plot of NIR light L 5 qq versus the 
binding mass M500, corresponding to the light and mass 
contained within the cluster radius r^QQ. We use Eqn 1 to 
estimate M 500 , with the cluster emission- weighted mean 
temperature Tx taken from the literature. The corre- 
sponding temperature scale is presented at the top of the 
figure. For our sample of 27 nearby clusters there is a 
reasonably tight correlation between cluster mass and K- 
band light 
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3 x 10 14 M Q 



0.69±0.09 



The best-fit results in this paper are obtained by mini- 
mizing the vertical distances in log space to the line from 
the data points, and no uncertainties in M500 are included 
in the fitting procedure. The la uncertainties are deter- 
mined by bootstrap resampling and refitting 10 3 times. As 
implied by the equation, the pivot point for the fit is at 
3 x 10 14 h^Q Mq, roughly corresponding to 4.8 keV. 

The best-fit correlation (solid line) differs from 
L500 oc M500 (dashed line) by ~ 3cr, indicating a 



mass-to-light ratio that increases with mass. In essence, 
low mass clusters produce relatively more light per unit 
binding mass than high mass clusters. Possible explana- 
tions include a higher star formation efficiency in low mass 
clusters (see §4.1) or a stripping process that liberates stars 
from galaxies more efficiently in high mass than in low 
mass clusters (e.g. Trentham & Mobasher 1998; Gregg & 
West 1998). As discussed in §1, we do not expect the 
if-band stellar mass-to-light ratio to be very sensitive to 
stellar population and star formation history variations, 
but gross differences in the stellar populations of low and 
high mass clusters could in principle contribute to the shal- 
low NIR light-mass relation found here. We return to this 
issue in §4. 

The fractional scatter in L500 about this best fit rela- 
tion is 28%, suggesting that If-band light does indeed 
provide a useful tracer of cluster mass. Possible sources 
of the scatter include: the variations in the cluster star 
formation and galaxy formation history, variations in the 
faint-end slope of the luminosity function a, the scatter 
(~ 17%) in the adopted M — Tx relation, the projection 
of unusually bright background galaxies onto the clusters, 
and the "deprojection" of the cluster luminosity (i.e. de- 
viations from spherical symmetry and from the assumed 
galaxy distribution model). In practice, the uncertainty in 
binding mass estimates will depend on the filter used to 
detect and characterize clusters. Crudely speaking, using 
the slope and amplitude of the scaling relation, for a given 
if-band light, one could predict the galaxy cluster mass 
with a statistical accuracy of ~ 45%. Of course this is 
likely a lower limit to the true mass uncertainties in a NIR 
cluster survey, because in practice one will not know the 
apparent virial radius #500 and cluster position a priori. 

Figure 5 contains the if -band mass-to-light ratio ver- 
sus binding mass Msoo- The trend of increasing if-band 
mass-to-light ratio T500 is best fitted by 



T500 = 47 ± 3 hro T Q 



5(10 



3 x 10 14 /ifo 1 Mi 



0.31±0.09 



(7) 



The Spearman correlation coefficient for the fit is 0.63, 
with a 1.6 x 10~ 3 probability of being consistent with 
no correlation. The average mass-to-light ratio over all 
clusters is 47 ± 3 hro Y © (short dashed line) and the aver- 
age for clusters with kTx > 3.7 keV (more massive than 
2 x 10 14 h^ M Q ) is 53 ± 3 h 70 T . A fit to the clusters 
with kTx > 3.7 keV produces a slope that is consistent 
with zero (0.03 ±0.12), which is one way of demonstrating 
that much of the leverage on the slope comes from the low 
mass systems. The 8 clusters that have kTx < 3.7 keV 
are of low redshift (0.016 < z < 0.057, z = 0.034). Be- 
cause visual inspections of these cluster fields do not show 
unusual behavior relative to higher Tx clusters (e.g. defi- 
ciency of galaxies), we believe the low T500 for these low 
Tx clusters is robust. 

In general, trends in blue/ visual band mass-to-light ra- 
tio with temperature have not been seen (e.g. Hradecky 
et al. 2000; David et al. 1995), but more recently both Bah- 
call & Comerford (2002) and Girardi et al. (2002) reported 
positive correlations between cluster blue band (V and B, 
respectively) mass-to-light ratio and temperature. Bah- 
call & Comerford (2002) observed T v oc T x 3±01 , which is 
somewhat shallower than the trend we see in our sample 
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(Tjc oc T£ 5±IU ); the best-fit of Girardi et al. (2002) sam- 
ple is T b oc M - 25 . Our results are statistically consistent 
with these other relations. Bahcall & Comerford (2002) in- 
terpret the dependence on temperature as an aging effect 
in the stellar population. If that is the case, then one would 
expect this trend to be less prominent in the NIR, where 
for old stellar populations the stellar mass-to-light ratio 
is relatively insensitive to population age. Because we see 
the trend in the -ftf-band, it is certainly possible that there 
is an underlying variation in star formation efficiency over 
cluster mass scales (we examine this possibility in §4.1). 

kT x (keV) 
2 4 6 8 10 




Fig. 5. — The cluster _R"-band mass— to-light ratio vs cluster mass. 
The solid line is the best-fit relation Eqn 7, the short-dashed line is 
the average of all clusters - 47 ± 3 hro^Qi and the long-dashed line 
is the average of hot (kT x > 3.7 keV, or M 50 > 2 X 10 14 /if 1 M Q ) 
clusters — 53 ± 3/i7oTq. The error bars mainly reflect the uncer- 
tainties in cluster binding mass estimated from Tx ■ 

Our results are consistent with the -KT-band 
mass-to-light ratio found in a detailed spectroscopic and 
2MASS photometric study of the Coma cluster (Tx = 
8.21 ± 0.16 keV): T = (68 ± 21) h 70 T Q (Rines et al. 2001, 
when converted from their projected mass-to-light ratio 
to the "deprojected" mass-to-light ratio adopted in our 
analysis). Interestingly, Kochanek et al. (2002) find 
a weakly decreasing mass-to-light ratio using 2MASS 
data. Evaluating cluster if-band light at r2oo, they find 



-200 oc 



Mr 



-0.10±0.09 



which is inconsistent with our re- 



sults at about the 4cr level. Their approach relies on a 
pseud-matched filter cluster finding algorithm using an 
NFW model with fixed virial radius and constant and 
a. The number of galaxies within the cluster virial region 
-/V*666 is estimated by rescaling the number observed using 
a fixed virial radius NFW filter. They show that iV*666 
is correlated with several estimators of binding mass. Al- 
though both our approach and theirs should, in principle, 
give consistent results, we believe our method offers several 
advantages for studies of the fC-band light - virial mass 
relationship. First, we start with a proven, low scatter 
X-ray estimator of cluster virial mass, and we use X-ray 
emission peaks as the cluster center. Thus, we know the 
size and location of our cluster virial region, whereas in 
their approach they have to pull these quantities from the 



2MASS data and bootstrap to an estimate of the total 
cluster light. Second, we use the NFW model only in 
deprojecting the cluster light, whereas their cluster de- 
tection algorithm is built upon the NFW model. Thus, 
our method would be less affected by deviations in some 
clusters from the NFW model. 

3.1. Cluster Mass-to- Light Ratio Constraint on Qm 

Cluster scale mass-to-light ratios T500, together with 
observations of the mean luminosity density of the uni- 
verse j, have often been used to estimate the mean mat- 
ter density ~p or matter density parameter S]m (e.g. Bah- 
call 1977a; Trimble 1987; Bahcall et al. 1995; Carlberg 
et al. 1996; Rines et al. 2001; Bahcall & Comerford 2002; 
Girardi et al. 2002, among others). Taking the mean 
mass-to-light ratio of the universe to be T un i V = ~p/j, 
then Qm ~ j^5ao/Pc if the cluster mass-to-light ratio 
is equal to the universal mass-to-light ratio. A potential 
weakness of this approach is that the galaxy populations 
inside and outside galaxy clusters differ significantly, with 
field galaxies being predominantly late type with a mix of 
young and old stellar populations and cluster galaxies be- 
ing predominantly early type with mostly old stellar pop- 
ulations (e.g. Dressier 1984). With these dramatic galaxy 
population differences, it is not clear that the star forma- 
tion efficiency within clusters would be representative of 
that in the universe as a whole. Nevertheless, an appli- 
cation of this technique using if-band light is particularly 
interesting, because it minimizes mass-to-light ratio dif- 
ferences due to the ages of stellar populations. Perhaps 
most interesting is that in combination with external con- 
straints on the matter density parameter, these results will 
provide insights into the star formation history in clusters 
as compared to the universe as a whole (see §3.3). 

Because the i-T-band cluster mass-to-light ratio is an in- 
creasing function of cluster mass, we will evaluate J7m by 
using both the mean T5oo, a /i = (47 ± 3) h 7 o Tq of all clus- 
ters, and the mean mass-to-light ratio of Tx > 3.7 keV 
clusters, Tsoojiot = (53 ± 3) h 7 Q T©. For the mean lumi- 
nosity density, we adopt the value obtained by Kochanek 
etal. (2001): J = {5.00 ± 0.52) x 10 s h 70 L G Mpc" 3 . Com- 
bining T5oo, a ;i or Tsoo./iot with j, with the critical density 
p c = 1.36 x 10 11 h 70 Mq Mpc -3 , we estimate the matter 
density parameter Q,m = 0.17 ± 0.02 (whole sample) or 
QjM — 0.19 ± 0.03 (hot cluster subsample). 

3.2. Total Baryon Fraction Constraint on VLm 

Another way of constraining the mean matter density of 
the universe from the clusters is by using the baryon frac- 
tion in galaxy clusters to estimate the universal baryon 
fraction. In combination with external constraints on the 
baryon density parameter, one can then estimate the mat- 
ter density parameter (e.g. White et al. 1993b; David et al. 
1995; Evrard 1997, MME). Hydrodynamical simulations 
indicate that the ICM is more extended in its distribu- 
tion than the cluster dark matter, and so corrections for 
this "depletion" effect must be included. In addition, hy- 
dro simulations indicate that clumping and substructure 
in the ICM enhance the X-ray emission relative to that ex- 
pected if the gas were smoothly distributed (MME; Math- 
iesen et al. 1999). A correction for this "clumping" effect 
must also be included. 
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Here we only consider the two dominant baryon reser- 
voirs in galaxy clusters: stars in the galaxies and the hot 
ICM (that is, we do not specifically include a contribution 
of the interstellar medium in galaxies, intergalactic stars 
or dark baryonic objects). We write the cluster baryon 
fraction as f b = {M star + M IC m )/M 500 , where M star and 
Micm are the total stellar and ICM mass inside r 50 o, re- 
spectively, and we have corrected the ICM contribution 
for depletion and clumping (see §2.3). We estimate M star 
using the if -band luminosity (see §2.2.3). If the corrected 
cluster baryon fraction then reflects the global one, we 
have Q M = 

kT x (keV) 
4 6 
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Fig. 6. — The baryon fraction as a function of cluster mass 
for the 13 2MASS clusters in the MME subsample. The long- 
dashed line is the mean of the 10 clusters that have kTx > 3.7 keV, 
fb hot = 0.160. The short-dashed line is the mean for all 13 clusters, 
J MH = 0.153. 

We show in Fig 6 the total baryon fraction in our clus- 
ter sample, which is an increasing function of cluster mass. 
The best-fit is 

„, f \ 0.148±0.04 

i«500 



f„ = 0.148^ 



0.005 



'© 



(8) 



for h 70 = 1. Note that this sample includes only 13 clus- 
ters - namely, those clusters in our ensemble that also lie 
in the MME X-ray flux limited sample of 45 clusters. The 
increase in the baryon fraction with cluster mass mainly 
reflects the dependence of the ICM mass fraction on clus- 
ter mass, because the ICM mass dominates over the stel- 
lar mass in galaxy clusters (we examine this issue in §4.2 
below). This correlation between ICM mass and temper- 
ature was noted by MME, whose study favors the trend 
M IC m/M 500 oc t£ 34±0 - 22 , which is a likely indication of 
preheating (e.g. Bialek et al. 2001) or perhaps variations in 
star formation efficiency (Bryan 2000). Because our deple- 
tion correction from hydrodynamical simulations is most 
appropriate for massive clusters, and because massive clus- 
ters are the systems that sample the largest portions of the 
universe, we take the baryon fraction in our massive sys- 
tems as most representative of the universal baryon frac- 
tion. In fact, the best-fit slope for the /& — Tx relation 
of the hot (kTx > 3.7 keV) clusters in this subsample is 



much shallower and consistent with zero: 0.054 ± 0.152, 
which may be an indication of the cluster baryon fraction 
asymptotically approaches the universal baryon fraction. 
Following the approach of the last section, we will present 
constraints on the density parameter that arise from the 
average characteristics of the entire subsample and the hot 
clusters in this subsample (10 clusters). 

The mean baryon fraction for the 10 clusters with 
kT x > 3.7 keV is /, 



b.hot 



0.160 ±0.006, and for the whole 

MME subsample we have J b all = 0.153 ± 0.006. We do 
not present the h scaling, because M star and Mjcm scale 
with the Hubble constant differently. With the value of 
fl b /j2 Q = 0.0457 ± 0.0018 from WMAP measurements, we 
estimate fijifj for the hot clusters in the subsample we 
get n M = 0.28 ± 0.03, while for the whole MME subsam- 
ple VLm = 0.30 ± 0.03. This value is consistent with the 
previous results using the same technique. Evrard (1997) 
deduced that Qj^ < 0.40 using the previous X-ray ob- 
servations. Based on ICM mass fraction measurements of 
45 k 30 clusters, MME and Ettori & Fabian (1999) con- 
cluded that < 0.32 & 0.34, respectively. Grego et al. 
(2001) estimated £Im ~ 0.25 from the SZE observations of 
18 clusters. Finally, with Chandra observations of 7 inter- 
mediate redshift clusters, Allen et al. (2002) deduced that 



M 



0.30 



-0.04 
0.03- 



3.3. Comparison of f2jf Constraints 

Our baryon fraction estimate of £Im is in excellent agree- 
ment with WMAP constraints from the CMB anisotropy: 
^M(h) = 0.28 ± 0.03 versus SI M (WMAP) = 0.27 ± 0.04 
(Bennett et al. 2003). Because the analysis that leads 
to the mass-temperature relation we use to estimate our 
cluster virial masses is so different from the X-ray imag- 
ing analysis required to estimate the ICM masses, the 
agreement between our baryon fraction estimate of S]« 
and the WMAP estimate suggests that our cluster virial 
masses and ICM masses must be fairly accurate. In 
that case we can use the baryon fraction Om estimate to 
constrain differences between the universal and the clus- 
ter mass-to -light ratio. Specifically, with the universal 
mass-to-light ratio T un i V — Slj^p c /j where SIm is given 
by cluster baryon fraction, we have T„„,„ = (78 ± 13) hm, 
whereas our hot cluster mass-to-light ratio is T500, hot — 
(53 ± 3) /170. Expressed slightly differently, our measure- 
ments imply that T U niv I ^boo, hot = 1-46 ± 0.21. It is well 
known that galaxy morphological mix and star formation 
rates are very different between the fields and the clusters 
(e.g. Dressier 1984; Lewis et al. 2002); our measurements 
thus provide an interesting indication of possible differ- 
ences in the star formation histories and star formation ef- 
ficiencies in these two dramatically different environments. 

4. STAR FORMATION AND ENRICHMENT 

Combined analyses of NIR and X-ray observations allow 
us to address several interesting questions concerning the 
thermodynamic history of baryons in the clusters: the star 
formation efficiency (§4.1), the mass fractions of stars and 
ICM (§4.2), and the metal enrichment of the ICM (§4.3). 
We discuss the implications of our results on it in §4.4. 

4.1. Variation of the Star Formation Efficiency 



Near-IR Properties of Galaxy Clusters 



9 





Tt nip 9 










Method 


Sample 




mass-to-light 


total (27) 


0.17 ±0.02 


ratio T 500 


massive (19) 


0.19 ±0.03 


baryon 


MME (13) 


0.30 ±0.03 


fraction ft> 


MME, massive (10) 


0.28 ±0.03 



Note. — Numbers in parentheses indicate the number of clusters in each sample/subsample. "massive" denotes clusters more massive than 
2 X 10 14 Mq. All quantities calculated assuming /170 = 1. 



Understanding the cosmic star formation history is of 
fundamental importance in understanding the formation 
and evolution of galaxies (e.g. Madau et al. 1998; Springel 
& Hernquist 2003) . Here we estimate the total stellar mass 
fraction for our sample of 27 clusters. This quantity can 
be regarded as a tracer of the total mass in stars formed in 
galaxies within the cluster halo over the course of its col- 
lapse. The total stellar mass for each cluster is converted 
from the total -ftT-band luminosity, using the derived mean 
stellar mass-to-light ratio (see §2.2.3 & Appendix). Be- 
cause NIR light is relatively insensitive to the star forma- 
tion history, our data should provide a robust estimate of 
the total stellar mass (within galaxies) in the clusters. 
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Fig. 7. — Fraction of halo mass that has been turned into stars 
by the present epoch in our cluster sample. The error budgets are 
dominated by the uncertainties in the cluster binding mass, which 
are due to the uncertainties in the X— ray temperature measurement 
and in the M — Tx relation (see Table 1). On the right axis is shown 
the cooled baryon fraction in the clusters (see §4.2 for definition). 
The dashed line indicates the universal cold baryon fraction (see 
§4.2). 

Figure 7 is a plot of the fraction of the halo mass which is 
in stars within galaxies for each of our clusters. This frac- 
tion varies from ~ 2.2% at low mass scales of lO 14 /^ 1 M & 
to ~ 1.2% at mass scales of 10 15 hj^ Mq. It is impor- 
tant to note that this fraction does not include any stars 
that have been stripped from their parent galaxies and are 
freely orbiting within the galaxy cluster potential, as has 



been observed in a some cases (Ferguson et al. 1998). Also 
shown in the figure is the best-fit relation 

Mstar - 1 R4+0.10 10 -2 ( M 500 \ (q] 

M 500 " 1 -°- 09 10 Ux 1014 M J (9) 
for hro = 1. Our observations suggest that the integrated 
star formation efficiency drops by a factor of ~ 1.8 over the 
mass range of galaxy clusters. Of course, it is also possible 
that the trend we observe arises because the reservoir of 
stars that lie outside galaxies is fractionally larger in more 
massive clusters (e.g. Trentham & Mobasher 1998; Gregg 
& West 1998). 

The decrease in star formation efficiency appears to 
be broadly consistent with theoretical expectations from 
state-of-the-art hydrodynamical simulations by Springel & 
Hernquist (2003). Their model includes radiative cooling 
and heating of the g multi-phase description of the 
interstellar medium, a self-regulating star formation mech- 
anism, and the feedback processes from supernova events 
and galactic outflows. They study star formation within 
haloes ranging from 1O 8 -1O 15 M0, and find star formation 
rates broadly consists with observations. In their models, 
the integrated star formation efficiency as a function of 
halo mass falls by a factor of 5 to 10 over cluster mass 
scales. The underlying cause is the less efficient formation 
of cooling flows in halos with virial temperatures above 
10 7 K (Springel & Hernquist 2003). However, what we 
observe is the star formation efficiency over all galaxy- 
mass halos that lie within our cluster halo at the present 
epoch. Because at early times the star formation efficiency 
in galaxies will be comparable for those galaxies in low 
mass and high mass cluster halos, we expect the theoret- 
ical prediction for M sta r / M^qo to suggest a much weaker 
trend. We are currently making a detailed comparison 
of our observational results with numerical star formation 
models. 

4.2. Variation of the ICM to Stellar Mass Ratio 

Another interesting quantity that has received much at- 
tention is the ICM to stellar mass ratio, Mjcm /M star , as 
a function of cluster mass. This quantity can be recast as 
the cluster cold baryon fraction, / CjC ; s = M star / (Mstar + 
Mjcm) = (1 + Micm /M s tar) , which is relevant in (nu- 
merical) studies of cosmic star formation (see Balogh et al. 
2001b, and references therein). In this section we examine 
these quantities for the MME subsample. 

Our joint analysis of NIR and X-ray observations shows 
that Mjcm /M star is an increasing function of cluster 
mass, varying from a factor of 5.9 for low mass clusters 
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to a factor of 10.4 for 
relation for /170 = 1 is 

Mjcm 



high mass clusters. The best-fit 



M. 



7 7+0.5 
'•'-0.4 



Mi 



-Am 



0.25±0.07 



^3xl0 14 M Q/ ' ! "" 
Over an order of magnitude in mass, there is roughly a two 
fold increase in Micm /M sta r) the existence of a trend in 
these data is significant at the 3.5 a level, providing con- 
firmation of the early work by David et al. (1990), who 
estimated the more massive clusters (Tx ~ 7 keV) have a 
~ 3 x larger Micm /M sta r ratio than do less massive ones 
(Tx ~ 3keV). Arnaud et al. (1992) found a correlation 
between Micm and stellar luminosities from early-type 
galaxies (E+SO); they estimated the mean ICM to stellar 
mass ratio as 3.6, which is much smaller than the values 
that we obtain. We note these two early studies eval- 
uated the ICM to stellar mass at different overdensities 
for different clusters. More recently, Roussel et al. (2000) 
found no trend when they examined Micm /M star at r2oo 
in their combined X-ray and optical sample of 33 clus- 
ters and groups. Interestingly, their analysis suggested a 
trend when measurements were made at T2oooj where less 
extrapolation of the X-ray data is required. 

x (keV) 
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Fig. 8. — The cold baryon fraction of the clusters. This is roughly 
inversely proportional to the ICM to stellar mass ratio. The solid 
line is the best-fit to the observed trend. The dashed line shows the 
universal cold baryon fraction. 

In Fig 8 we plot the cluster cold baryon fraction, the 
ratio of stellar mass to total baryon (star plus ICM) mass. 
As implied by Eqn 10, / CjC ; s is a decreasing function of 
cluster mass, indicating that the star formation efficiency 
is smaller in more massive clusters than that in low mass 
ones. This is in accordance with the conclusion of §4.1; 
in fact, from the ratio M s t a r/M^Qo one can infer the cold 
baryon fraction via / C)C / S = (M star /M 500 ) Cl M /Q b (Balogh 
et al. 2001b), which is shown in the right axis of Fig 7. Also 
shown in both Figures 7 & 8 (the horizontal dashed line) 
is the universal cold fraction f Cl univ = Qstar/^b, where 
Qstar = (2.7 ± 0.3) x 10 -3 /i^g 1 is obtained by an analysis 
of 2MASS data (Kochanek et al. 2001). We notice that 
in both figures the cluster cold fraction approaches the 
universal value with increasing halo mass. 



Our Figures 7 & 8 can be directly compared to the bot- 
tom and upper panel of Fig 1 in Balogh et al. (2001b), 
which was produced using previous observations of 48 
clusters and groups (Carlberg et al. 1996; Roussel et al. 
2000). In that study, it was argued that the first ap- 
proach (which utilizes the ratio Micm /M sta r) produces 
large scatter for f c , c ls, while the second one (that uses the 
ratio Mstar/Msoo) shows a much smaller scatter, and im- 
plies cluster cold fractions for all clusters considered are 
close to the universal value. However, measurements in 
our sample show this may not be the case. Our results 
may be more robust, because of the data homogeneity (all 
NIR photometry based on 2MASS), an improved treat- 
ment of stellar mass-to-light ratio, the use of observed 
X-ray M — T relation (as opposed to that based on nu- 
merical simulations), and the reliability of the ICM mass 
estimation (all X-ray data extend to r$oo). Studies of a 
larger sample will provide additional insights into these 
trends. 

4.3. ICM Enrichment 

Extensive investigations of ICM metal abundances have 
been carried out (see Arnaud et al. 1992; Renzini 1997 
for reviews), but until recently these studies relied on 
emission-weighted values and were therefore biased to- 
ward the abundance at the cluster center. Precision abun- 
dance measurement and radial abundance profiles in clus- 
ters have only been possible with imaging X-ray obser- 
vatories such as ASCA, BeppoSax, Chandra and XMM- 
Newton (e.g. Dupke & White 2000; Finoguenov et al. 2000, 
2001a; De Grandi & Molendi 2001, among others). 

In the bottom panel of Fig 9, we plot the iron yield, 
which is the ratio of iron (Fe) to stellar mass, as a func- 
tion of cluster mass in units of solar metalliticity: y c is,Fe = 

MFe/M star = (MiCM ZlCM,Fe + M sta r Z sta r,Fe) /M star 

(e.g. Arnaud et al. 1992). The ICM iron mass is calcu- 
lated by assuming the mean mass weighted iron abun- 
dance of the ICM to be Zicm,Fc — 0.21Z Q ^ e , which is de- 
duced from the metallicity profiles in cool core clusters and 
non-cool core clusters observed by BeppoSax (De Grandi 
& Molendi 2001, 2003 private communication). For the 
stellar iron abundance, we account for population trends 
in the galaxies much as we do in estimating the typical 
stellar mass-to-light ratio. We assume constant metal- 
licities Z star .Fe — l-2^o,_F e and 1.8ZQ^ e for ellipticals 
and spirals (J0rgensen 1999; Terlevich & Forbes 2002), re- 
spectively, and we weight them by the appropriate stellar 
mass-to-light ratio, luminosity functions of different mor- 
phological types, and the relative abundance of ellipticals 
and spirals. We detail our approach in the Appendix. The 
resultant Z star ,Fe is about 1.45-ZQ^e at Tx = 2 keV and 
about 1.27Z e>Fe at T x = 10 keV. ' 

The best-fit to the observed trend in our data is 

0.12±0.04 

j+0.10 
1.09 ' 



Vcls,Fe — 2.98_ ;q9 Z &t Fe 



Mi 



-.no 



(11) 



^3 x 10 14 M Q/ 

for /170 = 1. The iron yield is related to the iron 
mass-to-light ratio (IMLR, Renzini et al. 1993) which ap- 
pears frequently in recent literature (e.g. Finoguenov et al. 
2000): IMLR 

= Vcis,Fe^ star, where T s t ar is the mean stel- 
lar mass-to-light ratio. Previous studies using blue band 
optical photometry suggested that IMLR is generally in- 
dependent of cluster mass (Renzini 1997). The trend in 
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our sample for increasing y c is.Fe arises from the increasing 
ICM mass fraction, an improved understanding of the typ- 
ical ICM metallicity (De Grandi & Molendi 2001) and the 
decreasing stellar mass fraction. Fig 9 (lower panel) also 
shows that the iron yields in the cluster systems are very 
high relative to the approximately solar metal abundances 
observed in stars and the interstellar medium in galax- 
ies (Henry & Worthey 1999, and the references therein). 
It is unlikely that galaxies with a simple star formation 
history and typical initial mass function (hereafter IMF) 
could produce such high iron yields (although see Pipino 
et al. 2002). 
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Fig. 9. — The upper panel: the iron mass fraction in the clus- 
ters. The lower panel: the iron yield. The ICM iron abundance 
ZiCM,Fe = 0.21Z QiFe from De Grandi & Molendi (2001) is used, 
and the stellar contribution is calculated by the method stated in 
the Appendix. 

In Fig 9 (upper panel) we plot the iron mass fraction, 
f Fe = M Fe /M 500 , using Z Q . Fe = 1.814 x 10~ 3 (Anders 
& Grevesse 1989). Unlike the iron yield, which compares 
the iron mass to the stellar mass in each system, the iron 
mass fraction shows less apparent trend with respect to the 
cluster binding mass; the best-fit slope of the iron mass 
fraction is 0.048 ± 0.056, consistent with no trend at all. 
This nearly uniform enrichment of the baryon reservoir in 



clusters, combined with the high iron yield and the trend 
of decreasing star formation efficiency with cluster mass, 
may suggest that the primary enrichment mechanism in 
the ICM is not the stars that we observe in the cluster 
galaxies; instead, the intergalactic gas that becomes the 
ICM may be pre-enriched by some earlier stellar popu- 
lation that is at work before cluster galaxies form. This 
process appears to work with similar efficiency over the 
whole range of cluster masses. One possible mechanism 
for this uniform enrichment would be population III stars. 
An early stellar population of massive stars could con- 
tribute significantly to the metallicity - both because their 
likely top-heavy IMF and high metal yields (e.g. Larson 
1998; Heger & Woosley 2002) - and presumably enhance 
the entropy of the gas that eventually collapses to form 
galaxies and clusters. Another intriguing piece of cluster 
evidence for population III stars is that the observed ele- 
mental abundance ratios in the cluster ICM appear to be 
difficult to explain by reasonable combinations of Type I 
and Type II supernovae (Loewenstein 2001). 

4.4. Thermodynamic History of the ICM 

Clearly the intersection of cluster structure and the cos- 
mic star formation history is an intriguing avenue for fur- 
ther research. The star formation efficiency, the ratio of 
ICM mass to stellar mass and the total iron mass provide 
important insights into the thermodynamic history of the 
ICM. These results, taken together with existing obser- 
vations of other cluster scaling relations (e.g. Lx — Tx, 
Micm — Tx, Ri — Tx, where Ri is the X-ray isophotal size 
of clusters) provide an abundance of evidence that cluster 
ICM and stellar properties do not scale self similarly (e.g. 
MME, David et al. 1993; Mohr & Evrard 1997; Arnaud & 
Evrard 1999; Lloyd-Davies et al. 2000). One interpreta- 
tion has been that the gas distribution in low mass clus- 
ters has been altered relative to the expectations of the 
standard formation scenario by the addition of entropy or 
"preheating" by galaxy formation at early times (Kaiser 
1991; Evrard & Henry 1991; Cavaliere et al. 1997; Wu 
et al. 2000; Bialek et al. 2001; Muanwong et al. 2002). In 
this picture, prior to the cluster/group collapse, the inter- 
galactic gas is heated by, e.g., supernova-driven galactic 
outflows or AGN. The preheated, high-entropy gas thus 
forms a more extended distribution when it collapses into 
the cluster scale halos. The effect is less prominent in more 
massive systems, because in these massive systems the en- 
tropy increase in the ICM due to strong accretion shocks is 
much larger than that caused by preheating. This frame- 
work provides a means of explaining the observed slopes 
of the Mjcm — Tx, Lx — Tx and Ri — Tx relations, if 
preheating enhances the entropy by about ^100 keVcm 2 
(Bialek et al. 2001). This level of preheating appears to 
be consistent with the cluster "entropy floor" observed by 
Lloyd-Davies et al. (2000). 

Despite the successes of the preheating model in explain- 
ing observed cluster ICM properties, a generic problem 
with this model is the large energy required if carried out 
in the dense ICM in the centers of clusters or in the early 
universe (e.g. Bryan 2000; Wu et al. 2000). To circumvent 
these difficulties, Bryan (2000) proposed an alternative to 
the preheating model in which the star formation efficiency 
varies as a function of halo mass. The gas that cools to 
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form stars is the lowest entropy gas, and its removal from 
the cluster causes similar structural signatures to those 
produced in the preheating scenario (see also Voit et al. 
2002; Muanwong et al. 2002; Nath 2002). 

The NIR and X-ray properties of our cluster ensem- 
ble appear to provide problems for both the preheating 
and star formation efficiency scenarios. Namely, the direct 
evidence for a decreasing star formation efficiency with 
increasing cluster mass (Fig 7) lies outside the bounds 
of simple preheating scenarios, and the increasing total 
baryon fraction is at odds with simple models where the 
total baryon fraction is constant in clusters but the frac- 
tion cooling into stars is changing. Moreover, the high 
metallicity of the ICM together with the large ICM mass 
to stellar mass ratio strongly suggests that the source of 
enrichment (and so perhaps the entropy as well) requires 
something more than star formation with a typical ini- 
tial mass function. It appears likely that more complex 
models that include both varying star formation efficiency, 
preheating and enrichment at early times are necessary to 
explain these observations. 

5. SYSTEMATIC UNCERTAINTIES 

There are several key ingredients of our analysis: the 
cluster sample selection, the X-ray M — Tx relation, the 
ICM mass measurement, the total NIR luminosity, the 
mean NIR stellar mass-to-light ratio, and the ICM and 
stellar metallicities. These affect the derived NIR lumi- 
nosity, the stellar mass fraction, the cluster baryon frac- 
tion, and the iron mass fraction, as a function of cluster 
mass. In what follows we discuss the robustness of these 
ingredients. 

Cluster sample selection: Our cluster sample is taken 
partly from an X-ray flux limited sample and partly from a 
medley of clusters with available X-ray temperature mea- 
surements. Our results may depend on this selection, be- 
cause low mass clusters in our sample tend to be at lower 
rcdshift than the higher mass clusters. Any redshift re- 
lated systematic could then potentially masquerade as a 
mass trend. We examine this possibility in Fig 2, where 
we compare the luminosity function parameters for differ- 
ent mass and redshift ranges. As discussed in §2.2.2, after 
dividing the sample into two mass subsamples, we find no 
rcdshift trends in luminosity function parameters. 

M — Tx relation: The M — Tx determines the cluster 
mass and the virial radius. We adopt the observed rela- 



tion based on ASCA measurements of cluster temperature 
profiles (Finoguenov et al. 2001b). The self-similar cluster 
evolution scenario predicts the slope to be 1.5, which is 
consistent with the observed value we use. Using a slope 
of 1.5 (with a normalization chosen to be the same as the 
observed one at Tx = 6 keV) produces best-fit relations 
that are statistically consistent with the fiducial results. 
However, we notice that shifts in normalization can be 
important - in most of the cases crucial - to estimates of 
the cluster baryon fraction, mass-to-light ratio, and stel- 
lar mass fraction; however, normalization uncertainties do 
not affect the trends we see with cluster virial mass. 

ICM mass: All our ICM masses are from ROSAT PSPC 
X-ray surface brightness profiles that extend to or near 
rsoo, based on the adopted M — Tx relation. We also take 
into account the clumping in the ICM and depletion of the 
baryon fraction at r^QQ using results from hydrodynamical 
simulations (MME). 

NIR luminosity: We use the second release data from 
2MASS, which has accurate photometry and a small star- 
galaxy misidentification rate. Major factors that affect the 
estimates of the total NIR luminosity include: the cluster 
angular size (mainly by Tx & M — Tx relation), the faint- 
end slope a of the luminosity function, the contributions 
from the BCGs, and the method of background subtrac- 
tion. In the case that mass and light have the same dis- 
tribution, whether we choose rsoo or r2oo should not affect 
our results. Currently we are restricted by the sky cover- 
age of 2MASS second release data: we have NIR data only 
out to r 5 oo for some of the clusters, which increases the un- 
certainties in estimating light at r2oo- We have examined 
the light-mass relation at r 2 oo , and the best-fit slope is sta- 
tistically consistent with that at r 50 o ■ The faint-end slope 
of the luminosity function a has very little effect for those 
nearby clusters that 2MASS has probed deeply enough 
to measure a large fraction of light, but for clusters with 
larger redshifts, the effects are larger. Overall, though, the 
best-fit slopes calculated for several different values of a 
are consistent with the fiducial ones. Finally, the statis- 
tical and "annulus" methods for background subtraction 
give statistically consistent results. It is interesting to no- 
tice that, based on the observed if -band log TV — log 5 
relation, the probability that a foreground galaxy 2 mag- 
nitudes brighter than lies within the cluster region 
increases slightly with redshift for our sample; light from 
a chance superposition like this is subtracted statistically. 
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However, if there were any residual bias, the effect would 
be to add more light to higher redshift (which tend to be 
more massive) clusters, which would weaken the trend in 
mass-to-light ratio that we observe. 

Mean stellar mass-to-light ratio: We use this to esti- 
mate the stellar mass and the associated quantities. In- 
stead of using a varying mass-to-light ratio based on the 
(somewhat uncertain) spiral fraction and the (also un- 
certain) stellar mass-to-light ratio in elliptical and spiral 
galaxies, we could assume a constant mass-to-light ratio 
given by 2dF survey (Cole et al. 2001, see §2.2.3). Us- 
ing this constant mass-to-light ratio gives scaling relations 
that have slopes about 0.5 — la different from our fiducial 
ones. 

Mean ICM & stellar metallicities: These quantities are 
used to convert stellar and ICM mass into iron mass. The 
ICM iron abundance is measured to be relatively con- 
stant against Tx (De Grandi 2003, private communica- 
tion), and the value that we adopt is representative of 
the ICM as a whole. The stellar iron abundance, on the 
other hand, suffers difficulties that are similar to the mean 
stellar mass-to-light ratio, due to uncertainties in obser- 
vations and the relative spiral and elliptical abundance. 
If we assume solar abundance for stars in all galaxies, we 
obtain a yp e — M relation with slope ~ 1.5 a steeper than 
the fiducial one. 

In summary, as far as the best-fit results are concerned, 
all the effects considered above give consistent results with 
our fiducial results; the factors that give largest deviations 
are the faint-end slope of the luminosity function and the 
choice of cluster region. Other factors, such as the assumed 
cosmology in calculating the distance to the clusters, or 
the Malmquist bias (less than 3%), do not affect our re- 
sults. We therefore conclude that the results presented in 
the previous sections are robust. We look forward to bet- 
ter examining sample selection issues once the full 2MASS 
dataset is available. 

6. SUMMARY & PROSPECTS 

We have used the NIR data from the 2MASS survey to- 
gether with archival and published X-ray data to analyze 
the properties of the stellar and ICM baryon reservoirs in 
galaxy clusters. While X-ray temperature gives us reli- 
able cluster binding mass estimates, NIR light traces stel- 
lar population better than optical bands; the joint analysis 
makes it possible to examine whether stellar NIR luminos- 
ity can serve as a good cluster mass estimator. It also en- 
ables us to probe the star formation efficiency for clusters 
spanning an order of magnitude in mass. For a subsample 
of 13 clusters, we also measure ICM masses from X-ray 
data. We are therefore able to investigate the relative 
abundance of the main baryonic components in clusters, 
as well as the issues concerning the ICM iron enrichment 
processes. We review our conclusions here, and summarize 
the scaling relations in Table 3. 

First, the NIR properties of galaxy clusters appear to ex- 
hibit significant regularity. The if -band luminosity within 
r5oo is strongly correlated with the cluster binding mass, 
with an observed scatter of approximately 28% in if -band 
luminosity at a fixed cluster mass. Thus, the NIR clus- 
ter light is a good predictor of cluster mass. The if -band 
mass-to-light ratio increases from T 50 o ~ 33 /170 T Q for 
lO^/ifoM© clusters to T 500 ~ 6Sh 7a T Q for lO^/i^M© 



clusters. This trend must be accounted for when inter- 
preting NIR cluster surveys, where the NIR light is the 
primary indicator of cluster mass. However, we also no- 
tice that for clusters more massive than 2 x 10 14 /if 1 M© 
{Tx > 3.7 keV) there is little evidence for a changing 
mass-to-light ratio. 

Second, the cluster baryon fraction together with con- 
straints from the CMB observations on ill, provides a mea- 
sure of the matter density parameter Om = 0.28±0.03 (for 
^70 — 1) that is in excellent agreement with recent CMB 
anisotropy results. 

Third, the mass-to-light ratio of the most massive clus- 
ters can be used with a measure of the if-band luminos- 
ity density of the universe to estimate the matter den- 
sity parameter Q,m = 0.19 ± 0.03. Alternatively, adopt- 
ing the baryon fraction constraint on Qm, we show that 
the if -band mass-to-light ratio in the universe must be 
^univ — (78 ± 13)ft-7oT Q . This mass-to-light ratio is 
marginally consistent with that for our more massive clus- 
ters T500, hot = (53 ± 3)/i7oT Q , indicating that differences 
in the star formation history in these environments pro- 
duce present epoch mass-to-light ratios that differ by as 
much as ~30%. 

Fourth, the amount of if-band light per unit clus- 
ter binding mass is a factor of ~ 2 times higher for 
low mass clusters than for high mass clusters. Account- 
ing for a slight variation in the mean if -band stellar 
mass-to-light ratio for galaxies in low and high mass clus- 
ters, we find that the overall star formation efficiency de- 
creases with increasing cluster mass. This observation 
provides a new constraint on developing galaxy formation 
models. 

Fifth, the decreasing stellar mass fraction and increasing 
ICM mass fraction lead to a ratio of ICM mass to stellar 
mass that varies from 5.9 to 10.4 over the cluster mass 
range of 10 lA h^M Q to 10 15 h^M & . Thus, stars consti- 
tute roughly 14% of the baryons in low mass clusters and 
only 9% in high mass clusters. 

Sixth, the trends in stellar mass and ICM mass with 
cluster binding mass provide useful constraints on the ther- 
modynamic history of galaxy clusters. Specifically, pure 
preheating models are inconsistent with our observations, 
because preheating of the ICM would not introduce a 
decreasing star formation efficiency. In addition, simple 
cooling-oriented models that suggest a constant baryon 
fraction with a larger portion of the baryons going into 
stars in low mass clusters are inconsistent with our ob- 
servations, because the total baryon fraction is increas- 
ing with cluster mass in our sample. Models that include 
varying star formation efficiency and preheating or cool- 
ing to achieve a minimum entropy level in the ICM of 
~ lOOkcVcm 2 would likely reproduce the observed trends. 

Seventh, the iron yield per unit stellar mass is large 
and an increasing function of cluster binding mass. It is 
likely that unusually top-heavy stellar initial mass func- 
tions would be required to enrich the baryons to such a 
high level. Interestingly, the iron fraction (ratio of iron 
mass to cluster binding mass) is roughly constant. To- 
gether with the trend of decreasing stellar mass fraction, 
this constant iron mass fraction suggests that there may 
have been an early epoch of star formation and enrich- 
ment - perhaps before the epoch of galaxy formation - 
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that acted uniformly over the full range of cluster masses. 

Due to the incomplete sky coverage of the 2MASS sec- 
ond incremental release, our analysis is restricted to a 
small sample of 27 clusters; the shallow 2MASS photom- 
etry limits our study to low redshift clusters. We plan to 
expand the analysis presented here to a much larger sam- 
ple of local clusters, as well as to clusters at higher redshift 
using deeper photometry. Together, these studies will lead 
us toward a more complete understanding of the evolution 
of galaxy clusters and their baryon reservoirs. This under- 
standing will undoubtedly contribute to interesting cosmo- 
logical studies and an improved understanding of the star 
formation history of the universe. 
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APPENDIX 

Here we describe our estimate of the mean stellar mass-to-light ratio for each cluster. As a simplified model, we divide 
galaxies into early and late types and calculate the average stellar mass-to-light ratio in the two types. We employ 
published estimates of the stellar mass-to-light ratio as a function of mass and luminosity in both ellipticals and spirals. 
With published estimates of the spiral fraction as a function of cluster mass and of the luminosity function for ellipticals 
and spirals, we then calculate the luminosity-weighted mean stellar mass-to-light ratio. We describe each of these steps 
in turn. 

For the if-band stellar mass-to-light ratio in ellipticals, we use the results of a dynamical analysis of 21 luminous ellip- 
ticals (Gerhard et al. 2001). Specifically, from their Fig 13 we estimate the central mass-to-light ratio T e (L) as a function 
of galaxy luminosity, using the typical color for ellipticals B — K w 4.1 (Pahre 1999). The central mass-to-light ratio is a 
good estimate of the stellar mass-to-light ratio, because the central regions of ellipticals are dominated by stellar mass. 
For the if-band stellar mass-to-light ratio of spiral galaxies, we use the results of Bell & de Jong (2001, see their Fig 1), 
who construct models that describe many characteristics of spirals in the Ursa Major cluster. 

We adopt the K-hand luminosity functions (characterized by the Schechter parameters a and £*) from a K-h&nA 
2MASS study (Kochanek et al. 2001). With subscripts e & s denoting "elliptical' & "spiral", respectively, the mean 
stellar mass-to-light ratio is 

n e L*, e T e (y) e'Vy^dy + n s L*, s f™^ T s (y) e-Vy°*+Hy 

Tstar ~ n - L *<- -d.e e-yy^+idy + n s L^ s e-Vy^dy ' 

where yi ow ,i = Li ow /L*^ (i — e, s), Li ow corresponds to M\ ow = —20, the luminosity cutoff (see §2.2.3), the galaxy number 
densities n e & n s are determined from the total number of galaxies observed n' tot = , L 0* (L/L*) a cxp(—L/L*)d(L/L*) 
and the spiral galaxy number fraction f s : n e = n' tot (l — /„)/ e~ y y a ° dy , n s — n' tot f s / e~ v y as dy. As explained 

earlier, the integration limits are chosen to avoid integrating over the faint-end of the luminosity functions whose shape 
is unknown and where the behavior of galaxy T(L) is also unknown. 

We use early published estimates of the spiral fraction to measure the trend in f s with cluster Tx (Bahcall 1977b; 
Dressier 1980; Dressier & Shcctman 1988). For the 28 clusters from these data sets with published Tx, we regard all 
galaxies designated as cD, E and SO as ellipticals, and we ignore irregulars. The resulting spiral fraction decreases with 
cluster temperature. The best-fit to the f s — Tx relation is log/ s = —0.5 log(Tx/4kcV) — 0.4. For each cluster in our 
sample we calculate f s (Tx), then use the above expression to calculate the mean stellar mass-to-light ratio. The resulting 
T" star varies weakly with Tx (from ~ 0.7T Q at 2keV to ~ 0.8T Q at lOkeV). Curiously, the value at low mass end is very 
close to that obtained by Cole et al. (2001), using a Kennicutt (1983) IMF (see Fig 3). 

Let us point out some caveats for the above approach. First of all, the galaxy mass-to-light ratios are not solely 
from observations of cluster galaxies (elliptical mass-to-light ratio are observed but not for galaxies in clusters; spiral 
mass-to-light ratio are from models built from observations of cluster galaxies, see above), and the luminosity functions 
used to calculate the mean mass-to-light ratio are based on all the galaxies with cz > 2500 km s _1 , which includes both 
cluster and field galaxies (Kochanek et al. 2001). We simply assume that these "field" luminosity functions are applicable 
to cluster environments. Finally, the spiral abundances presented in Bahcall (1977b) and Dressier (1980) are based on 
visual/blue band observations (which may be larger than the values obtained when observed at redder bands), and are 
not restricted to a fixed fraction of virial radius or a fraction of the luminosity function comparable to that we use in our 
2MASS study; thus, the uncertainty in f s is therefore probably significant and should be examined directly in the 2MASS 
data. 

Following a similar line of reasoning, we calculate the mean stellar metallicity based on the observed iron abundance 
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for stars in elliptical and spiral galaxies Z e & Z s : 

_ n e £., e Z e (y) T e (y) e^y^dy + n s L., a Q„„ Z s (y) T s (y) e^y^dy 

Zstar - FC ~ n e L^ e C T e (y)e- y y a ° +1 dy + n s L.. s f°° ' T s (y) e-vy^dy " 

For stellar iron abundance in ellipticals we use the results from a study of 115 E/SO galaxies in Coma (J0rgensen 1999). 
This analysis showed that there is no clear correlation between [Fe/H] and galaxy mass or luminosity. The median value 
of [Fe/H] is ~ 0.1, which corresponds to Z e ~ 1.2Z Q ^ e . For spiral metallicity, we examine measurements for 14 spirals 
in cluster/group environments (Terlevich & Forbes 2002). [Fe/H] of these spirals does not show a clear correlation with 
galaxy luminosity; the mean [Fe/H] ~ 0.25, or Z s ~ 1.8Z Q j? e . Assuming these metallicities are representative for all 
spirals and ellipticals, the resultant mean metallicity Z sta r Fe is about 1.45 at Tx = 2 keV and about 1.27 at Tx = 10 
keV. 
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